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The generating function for all antisymmetric characters of a Wilson loop matrix in SU (N) Yang 



O | Mills theory is the partition function of a fermion living on the curve describing the loop. This 



generalizes to fermion subsystems living on higher dimensional submanifolds, for example, sur- 
faces. This write-up also contains some extra background, in response to some questions raised 
during the oral presentation. 
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1. Introduction 

In D Euclidean dimensions, at N = °° SU (N) pure YM simplifies [1], but this has not yet been 
turned into a quantitative tool. We are trying to improve on this situation. 

Let & be an observable characterized by one single scale / in 4D. For IA^ <C 1, {@)n=°o can 
be computed by summing up contributions from planar Feynman diagrams. Moreover, there are 
good reasons to believe that the series in log ^ Ajv ^ converges [2]. Am is the analogue of Aqcd for 
pure SU (N) gauge theory. For IA^ 3> 1, for a class of specific ^'s, {&)n=oo can be expanded in 
j£- using an effective theory based on free strings [3]. 

How are the two regimes connected ? 

Some £?'s have a narrow crossover as / changes from I An < 1 to IA^ > 1 at finite N, becoming 
a "phase transition" atN = <*>. Examples of such will be presented later - in these examples the 
operators turn out to be non-local [4]. Such large N phase transitions tend to fall in Random Matrix 
universality classes. The hope is to exploit this to connect field theoretical perturbation theory to 
effective string theory. 

2. Loops 

2.1 Wilson loop operator - ignoring renormalization 

Let R be an irreducible SU (N) representation of dimension oIr. We define the nonlocal matrix 
valued operator: 

&>e i ** AR - dx = a R (tf) (2.1) 

& denotes ordering round the closed curve . The starting/end point of the ordering is suppressed 
because we shall build out of Q. R (^) other operators which will all be independent of the choice of 
this point on ^ . The Wilson loops are the following averages: 

W R (tf)=ti(n R (tf))/d R (2.2) 

Typically, one focuses only on the fundamental representation, which we denote by R = 1 or 
R = f. We depart from this by considering the collection of all fully antisymmetric representations, 
R = k,N-k;k=l,..., [fj.The generating function for the W k (V) is (det(z+Q/(«0)) = g(z,?f), 
a palindromic polynomial of rank N in z. This permits us to introduce a set of quantities we loosely 
refer to as the "eigenvalues" of ni(^), the roots of Q(z,^). As the loop is dilated, its minimal 
area A grows and one can extract ^-string tensions from W^ffi) ~ exp(— GjA). 

The following scenario is plausible: A large Af "phase transition" (meaning: a non-smooth 
change in the dependence on z) separates small from large loops and occurs in D=2,3,4. As N — > °° 
a continuum density of "eigenvalues", roots of Q(z, e ^'), develops, supported on \z\ = 1, gapped at 
z = 1 for small loops, and uniform for infinite loops. Hence, for very large loops the deviation of 
the density from uniformity is controlled by [N/2] exponents e~° kA , dominated by e~ c f A . Here we 
imposed an assumed parity invariance. 

Close to critical loop-size, and for z close to 1, there is a universal description common to all 
dimensions. The case D = 2 is exactly soluble so the universal form is known [5]. This phase 
transition is seen only in Q{z, c <o), but not in the individual WkS and is therefore different from the 
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familiar Gross-Witten phase transition [6] . The universality provides an economic parametrization 
of the short-scale to long-scale crossover in Q{z, c tf) for 1 <C N < °°. 

For this to be meaningful we need to renormalize Q(z, c &)- After this is done, we can use lattice 
gauge theory simulations to validate the scenario as a property of continuum large N SU (N) gauge 
theory. 

2.2 Renormalization of Wilson loops 

Our bare observable is Q(z, c rf) = (det(l + zQAff))) . The following representation is an exact 
identity for an Clf given by a fixed N x N unitary matrix of unit determinant [7]: 



Here, z = e~^ 1 ', o parametrizes *«f by x(o) with [d a x^(a)] 2 = 1 and / is the length of *£. The 
Grassmann variables vK<t), Y( a ) obey anti-periodic boundary conditions. The gauge potential 
along the curve is projected from the D-dimensional gauge potential by a(a) = (jc(ct)) . 
We rewrite out bare observable accordingly: 



Now we employ power counting to determine the needed counterterms (ct-s) for renormaliza- 
tion [8]. The parameterization choice of the curve makes the curve parameter of dimension inverse 
mass, [a] = — 1, which in turn makes the Grassmann fields unitless, [Xfr, yf] = 0. Terms containing 
more than one xj/y pair and a single derivative can be eliminated by a field redefinition. This leaves 
the non-redundant ct-s, [xj/Xj/f , k = l,..N. The symmetry \jf — > y/, \]r — > — > A* can be 



This leaves exactly the amount of freedom needed to arrange for the inequality ^Wr(^) < 1 to 
hold for all antisymmetric R, after renormalization. So long as this is true, some further adjustments 
(not finely tuned) of the finite parts of the coefficients of the ct-s will render Q(z, 'if) palindromic 
with all roots on \z\ = 1. To be sure, in 4D and 3D, the short distance singularity that one expects 
from one gluon exchange even after the removal of the perimeter divergence, is eliminated by the 
integration round the contour and one can set the expectation values to their classical values for 
zero size loops. In 4D for example, the one gluon exchange answer could be regularized to give a 
nontrivial conformally invariant functional of the loop shape [9], but we do not do this here. 

The divergences one has to eliminate are linear in 4D and logarithmic in 3D. There are no 
divergences in 2D. All our loops are assumed to be differentiable, so the corner singularities of the 
4D case do not occur. 

2.3 A Dirac operator associated with the large jV transition 

The spectrum of D\ = d a — ia(o) will have a gap for small loops and will be gap-less for 
large loops. There is an analogy to spontaneous chiral symmetry breaking and its connection to 
chiral random matrix theory. This analogy is exploited in our generalization to surfaces. 




(2.3) 




(2.4) 
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2.4 More about the divergences 

The [j] ct-s are necessary to eliminate the [j] perimeter divergences associated with the 
distinct A^-ality representations, not counting conjugate ones. Physically, the ct-s represent the 
arbitrary amounts of thickening the distinct ^-strings need. 

On the lattice one "thickens" the gauge fields the fermions see. For example, in APE smearing 
[10] one defines a recursion starting from the link matrices Uu(x) at step n = and proceeds from 
step n to step n + 1 by: 

X^ +l \x;f) = (l-\f\)U^\x;f) + ^-^ ■ 



2{d-\) WW) 

, , = 1 : =— (2-5) 
[X^\x;f)n^\x;f) 

in) 

£ („), . denotes the "staple" associated with the link U u (x;f) in terms of the entire set of 

Uy 1 \y;f) matrices. Writing U^\x;f) = exp(iA^ (x;f)), and expanding in one finds [11], 
in lattice Fourier space: 

4" +1 W)=Lv(?)4"W) (2.6) 



with 



where = 2sin(^-) and 



V(ff)=/(«)(*Mv-^) + ? ^ (2-7) 
q z q L 

f^= l -2(h-/ (28) 

The iteration is solved by replacing f(q) by f"(q), where, for small enough /, 

/"(<?) ~e~5(^ 2 (2.9) 

One can think about f"(q) as a gauge invariant, Euclidean, form factor in a particular gauge; this 
form-factor is attached to the gauge field - fermion pair vertex. 

One can take a trivial continuum limit [12] on the step number and fn is then seen to play 
the role of the fifth "time", %, in a Langevin equation where the noise has been set to zero. The 
appropriate D+ 1 dimensional gauge covariant set of partial differential equations is 

F 5v =D^ v , (2.10) 

where we use continuum notation also for the gauge fields and have set % = x$. These equations 
have only an 0(D) invariance, not 0(D + 1). They are first order in ^5 and second order in , jJ. = 
1,2,3,4. Gauge fixing to A5 = and going back to t, we get 

dA tl 

DnFuv (2.11) 



dr 



with a residual 4D gauge invariance - see [12] for references to earlier work on similar equations. 

The parameter T tells us, in units of length-squared (2.9), the rough amount of regularizing 
"fat" that has been added to the loop. All ultraviolet divergences of Wilson loop operators are 
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eliminated as soon as x > 0. These would include corner divergences if there were any; on the 
lattice there always will be. 

We take x proportional to the physical perimeter of the loop squared for small loops, and 
approaching a constant length squared (relatively small relative to ) for large loops. This is 
done in order to remain sensitive to the short distance structure of the theory for small loops while 
allowing large loops to have only a modest amount of thickness. The confinement of the fat loops 
will still be governed by exponential terms with the same Ok one would have extracted from bare 
Wilson loops, which were not fattened. 

3. Surfaces 

The need to smear the loops and to adjust a number of order N of finite parts renders the 
renormalized Q(z, c <o) observable somewhat inelegant. This leads us to the idea of looking at less 
singular objects than lines - less singular means that more displacement directions are allowed. 
The simplest step in this direction is to replace the curve by a 2-dimensional surface E, described 
by jc^(cj) (da, a = 1,2). We put massive Dirac fermions on E, which is characterized by a single 
scale I. At zero mass, one has chiral symmetry and the fermionic determinant is the exponent of the 

dx 

Polyakov-Wiegmann action [13]. The gauge connection on £ is a a = A^(x(o))-^- . The massless 
Dirac operator is £>2(£) = Ya[d<j a — 'fla(d)] and our new observable becomes 



The possible counterterms are built out of currents on the surface: J J a (o) = \j/(a)YaT J Y(<T), 
Ja = Y(o)YaW(o)- So, only two ct-s, Jzfi = J J a J J a , Jz?2 = JaJa are now available. In 4D one 
has logarithmic divergences, and Jz?i will be generated, but in 3D no ct-s are generated because the 
theory is superrenormalizable. 

At N = °° a S^SB transition is possible, depending on I; for a finite surface, this is a large 
N phase transition, analogue to the one occurring in the smeared Wilson loops. For finite ,/V this 
transition disappears. 

The lattice version of <2(/i,E) is straightforward to construct: we used the overlap operator for 
D2 [14] to ensure a precise definition of the location of spontaneous chiral symmetry breaking on 
the lattice. 

4. Results in 3D 

As said, no extra renormalization necessary in 3D. We set pL = and N = °°. For an infinite 
plane we find {xj/y}^*, = 0.29(l)y / ay, where ay is the k = 1 string tension in the pure YM theory. 
For a cylinder with square base of side s we establish that spontaneous chiral symmetry breaking 
occurs when s > s c , while for s < s c chiral symmetry is preserved. For s > s c with (s — s c )/s c <C 1 
we have {yy)n=^{s) x (s — s c ) l l 2 . The critical size is found to be s c = 1.4(l)/y / a/. 





(3.1) 
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5. Outlook 

In 3D the definition of the curve related observable that exhibited the large N phase transition 
required extra regularization which we implemented as smearing. Moving from curves to surfaces 
eliminated this smearing. In 4D there are further obstacles which we have to learn how to overcome. 
However, the number of ct-s is independent of N now. 

In 3D, we described above two unambiguously defined N = oo scales, (y/y/) and s c . We are 
looking for analytical ways to estimate these scales. 
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